M = 40; N = 1e5;
x_n = randn(1, N);
x_n = 2*x_n + 16;
figure(1);
hist(x_n, 40);
hold on;